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SPACE VLBI FOR GEODESY AND GEODYNAMICS 


1.1 Introduction 

In the 1990’s dedicated radio telescopes will be launched into Earth orbit and will be 
integrated in the ground-based VLBI networks. A straightforward extension from present 
ground-based VLBI to space is called space VLBI , which uses radio-antennas in space. 

The first space VLBI mission that will be implemented in the near future is 
RADIOASTRON. Its trajectory is planned to be 3,000 km in perigee and 69,000 km in 
apogee. The satellite will carry a 10 m antenna. It will be launched in 1993. The second 
project is a Japanese orbiting VLBI mission called VSOP. Its trajectory is 1,000 km in 
perigee and 20,000 km in apogee. The expected launch date is 1995. Both projects are 
now in progress. Most probably a combined use of RADIOASTRON and VSOP 
observations will be performed in the future. 

Space VLBI observables may be useful to improve Earth’s gravity field and in the 
unification and connection of ref erence frames inherent in the space VLBI technique. The 
space VLBI missions offer and will provide new types of satellite observables (VLBI time 
delay, delay rate and differential VLBI tracking data) with high accuracy for these potential 
applications. 

In the simplest version of the space VLBI technique, one station in orbit observes in 
conjunction with a second station on the ground. However, in practice we have a number 
of networks of ground antennas observing the common celestial radio sources 
simultaneously with a conventional VLBI technique. Moreover, joint observations of two 
or more space VLBI satellites will supposedly be performed in the future'. Therefore, in 
our research work a combined use of simultaneous space and ground-based VLBI 
observations are considered from the geodetic and geodynamic points of view. 

1.2 Purpose of the Research 

The space VLBI observables may be useful to improve the Earth’s gravity field and in 
the unification and connection of reference frames inherent in the space VLBI technique. 
The space VLBI missions offer and provide new types of satellite observables (VLBI time 
delay, delay rate and differential VLBI tracking data) with high accuracy for these potential 
applications. Therefore, we started to investigate potential possibilities of space VLBI 
missions for these areas. The purpose of the corresponding research work is to explore the 
feasibility of these potential applications and to provide sufficient background for the 
inclusion of space VLBI observables in geodetic data processing programs (e.g., 
GEODYN). The final goal of the research from the point of view of geodetic science 
should be to figure out whether or not space VLBI can provide real improvements in the 
two potential applications in geodesy and geodynamics. 

The purpose of this research work is, as a first step, to derive the mathematical models 
of space VLBI observables suitable for least squares covariance analysis as well as to 
explore estimability problems inherent in the space VLBI system, including a detailed rank 
defect analysis and sensitivity analysis. An important aim is to carry out a comparative 
analysis of the mathematical models of the ground-based VLBI and space VLBI 
observables in order to describe the background in detail. Developing computer programs 
in order to check the relations, and for error assessment and sensitivity analysis later on, is 
also an important purpose of this work. 
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1.3 Estimabillty of Geodetic and Geodynamlc Parameters from Space VLBI 

In order to investigate the estimability of different geodetic and geodynamic parameters 
from the space VLBI observables, the mathematical models for time delay and time delay 
rate observables of space VLBI have been analytically derived along with the partial 
derivatives with respect to the parameters. Rank defect analysis has been carried out both 
by analytical and numerical testing of linear dependencies between the columns of the 
normal matrix thus formed. Definite conclusions have been formed about the rank defects 
in the system. 

1.3.1 Mathematical Models. The mathematical models for ground-to-space VLBI 
observables have been defined through equations formed by suitable modifications to the 
equations relating the ground-to-ground VLBI observables to the parameters. Our main 
purpose is not to derive explicit observations equations, but rather to develop qualitative 
expressions to demonstrate the relationship between the observables (time delay and delay 
rate) and the parameters (station and source coordinates, orbital parameters, earth rotation 
parameters — ERP’s, and clock parameters), and to explore the estimability problems 
inherent in the system. 

The mathematical models for time delay and delay rate observables for ground-to-space 
and for space-to-space VLBI observations have been derived and are included in a detailed 
report by Jozsef Adam (Estimability of Geodetic Parameters from Space VLBI 
Observables, Department of Geodetic Science and Surveying Rep. 406, in press). 

1.3.2 Rank Defect Analysis. To carry out the analysis of rank deficiency of the 
normal matrix for estimating the set of parameters from the time delay ground-to-space 
observable, partial derivatives of the observables with respect to the parameters are derived 
analytically as a first step. Examining the analytical expressions for the partial derivatives, 
the four linear dependency relations between the columns of the design matrix have been 
derived as given in the abovementioned report. 

Thus from the analytical method, the rank deficiency of the system is found to be four. 
In order to confirm this, numerical testing was carried out using software developed/ 
modified/obtained from various sources. The stations and radio sources and their 
coordinates, the orbital parameters of the space VLBI satellite and the basic configuration 
used in the test computations are tabulated in Tables 1-4 of the report by Adam. 

To determine the rank deficiency of the normal matrix, eigenvalue analysis is done, the 
results of which for both the short-arc and full orbit cases are given in Table 5 of the report 
by Adam. These are also confirmed by doing Cholesky lower-upper triangular matrix 
decomposition of the matrix. It is obvious from the results that the rank deficiency of four, 
as derived analytically, is confirmed from these numerical tests. 

1.3.3 Conclusions. From the analytical derivations and also from the numerical 
tests carried out, it is concluded that the rank deficiency of the normal equation systems in 
the case under investigation is four. A geometric interpretation for this rank deficiency can 
be given as follows: A linear dependence between various partial derivatives of the ground 
station coordinates and the earth rotation parameters indicate a rank deficiency of three due 
to a lack of absolute orientation of the network of ground stations with respect to the true- 
of-date celestial frame which cannot be sensed by the observables. The fourth rank 
deficiency, a linear dependency between partial derivatives with respect to the right 
ascension of a radio source and the ascending node of satellite orbit, is due to a lack of 
reference direction for the true-of-date celestial frame, i.e., due to the observables being 


2 


insensitive to the orientation of the true-of-date inertial frame in right ascension. The origin 
of the reference frames is fixed by the given Keplerian orbital elements at some reference 
epoch To, thus rendering a rank deficiency of only four. For more detail, please refer to 
the report by Adam. 

1.4 Sensitivity Analysis 

After concluding the rank defect analysis of the system, numerical testing was taken up 
to determine the sensitivity of the space VLBI time delay observables to the determination 
of-different parameters of geodetic and geodynamic interest, including the ERP’s, orbital 
parameters, etc. Computer programs have been developed/modified/obtained and used for 
such numerical testing. 

As a first step, Cartesian coordinates of the space VLBI satellite at specified epochs and 
coordinates of the quasars are obtained and converted to Keplerian orbital elements. Input 
files with UT and CAST epochs at a given interval are created, and ERP’s are predicted 
using the IERS formula for the prediction. The partial derivatives are numerically evaluated 
using a computey)rogram developed to obtain elements of design matrix A from which the 
normal matrix A* A is computed, assuming unit weights. 

Different combinations of number of epochs, number of stations and number of radio 
sources are used to investigate the effect of these factors on the sensitivity. These 
computations have been carried out for both the proposed space VLBI systems: VSOP and 
RADIOASTRON. Different combinations of parameters to be estimated have also been 
used in order to study the sensitivity to different combinations of parameters. 

The numerical results of the different tests carried out indicate that the accuracy of the 
ERP’s determined from the minimum configuration of three observing stations, three 
quasar sources and ten observations epochs, for an observation period of about six hours is 
as low as a few tens of milliarcseconds, which improves to the order of about 0.3 
milliseconds for a configuration of 15 stations, nine sources and 2880 observations 
epochs, with an observation period of about two days. 

It must be noted that these accuracy estimates are most optimistic, assuming a simple 
Keplerian orbit, no errors in the parameters held fixed, and before fixing the rank defect. 
These assumptions are only gross approximations and cannot be considered realistic. A 
more detailed sensitivity analysis, by fixing rank defects analytically and by fixing 
parameters and including error propagation and also by removing some of the 
approximations mentioned above, is in hand. The results from the preliminary simple 
analysis indicate accuracies marginally better than those obtained by the existing techniques 
for the determination of the ERP’s; however, a more detailed investigation may be required 
before reaching definite conclusions. 

1.5 Error Analysis for Estimation of Earth Rotation Parameters 

From the sensitivity analysis carried out to analyze the effect of the different 
parameters, e.g., the number of observing stations, number of radio sources, number of 
observation epochs, total observation time, etc. on the accuracy of earth rotation parameters 
to be estimated from the space VLBI observables, the optimum number of these parameters 
were found to be: 

Number of observation epochs 1800 

Number of observing stations 15 

Number of radio sources 3 

Total observation time 45 hours 


3 



With these optimum number of parameters, a detailed error analysis to estimate the 
effect of errors in the parameters on die errors in the ERPs estimated from this data, using 
the simplified mathematical model, was taken up. A priori variances were assigned to the 
parameters: station coordinates, radio source coordinates, orbital parameters (Keplerian 
elements of the satellite orbit), and the clock parameters, and these errors were apropagated 
into the ERPs using the simplified mathematical model. 

The necessary relationships and models for error propagation were derived and 
computer programs developed for computing these effects numerically, using the simulated 
orbit, the necessary data and the assumed a priori variances. Different values were chosen 
for the a priori variances assigned to different parameters, from the most optimistic case to 
the worst case. 

The results obtained through this detailed error analysis were tabulated and plotted and 
were presented in the poster papaer titled “Space VLBI Error Analysis Results for Earth 
Rotation Parameters,” presented at the NASA CDP Investigators’ Meeting, Greenbelt, 
Maryland, October 25-26, 1990. (See the Appendix.) 

From the results obtained, it can be seen that in the most realistic case — Case B (Page 
17 of Appendix A) — the error in the ERPs estimated from the space VLBI observables is 
about 0.7 mas in ERPs X p and Yp, and about 0.06 ms in the change in the earth’s rotation 
rate (UT1-UTC), which are still marginally better than the accuracy obtained by present 
techniques. Hence it is considered worthwhile to continue these investigations and to take 
up a detailed simulation analysis for the feasibility study of determination of the ERPs from 
the space VLBI observables. It can also be seen from the results of analysis of effect of 
orbit errors on ERPs (page 15 of Appendix A) that errors in orbit, holing the station 
positions fixed, are not propagated effectively to the errors in ERPs through the present 
model. Hence a detailed investigation of this also is in hand. 
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